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TECHNICAL NOTE 2698 


THEORETICAL ANALYSIS OF HYDRODYNAMIC IMPACT OF 
A PRISMATIC FLOAT HAVING FREEDOM IN TRIM 
By Robert W. Miller 


SUMMARY 


Equations which include freedom in trim are derived for Hydrodynamic 
impact of a non- chine -immersed, prismatic float forebody having a V-bottom 
and a transverse step. These equations are an extension of previously 
published fixed-trim theory, and a method of solution is indicated by 
which time histories of vertiqal, horizontal, and angular displacement, 
velocity, and acceleration can be obtained. 

Comparisons of specific solutions of the equations with corresponding 
fixed- trim solutions are presented. The trends and deviations noted are 
similar to those exhibited by a like comparison of experimental data for . 
free- and fixed- trim impacts. 


INTRODUCTION 


Previously published hydrodynamic impact theory (references 1 to 3 ) 
has been based on the concept that the flow about an i m mersing seaplane 
float or hull is a two-dimensional phenomenon occurring in transverse 
planes fixed in space and oriented normal to the keel. The total force 
on the float is obtained by summing the reactions in the individual flow 
planes in contact with the float and applying an aspect-ratio correction 
factor to account for the end-flow losses which exist in three-di m ensional 
flow. This theory has made use of the si nullifying assumption that the 
trim remains fixed throughout the impact. Experimental checks of this 
fixed-trim theory for both model and full-scale hulls have been presented 
in numerous reports and some evaluation of the empirical factors involved 
has been conducted. 

The present investigation was initiated in order to obtain a method 
for determining the effect of freedom in trim on loads, moments, and 
motions during hydrodynamic impacts of a non-chine-immersed, prismatic 
float ' forebody having a V-bottom and a transverse step. As in the pre- 
vious theoretical presentations the wing lift was assumed to be equal 
to the model weight and the aerodynamic moments were neglected in order 



2 


NACA TN 2698 


to simplify the problem, although they might have had an appreciable 
stabilizing effect. 

The present paper is primarily concerned with the derivation of 
equations describing the vertical, longitudinal, and angular motions and 
accelerations of the previously described float during hydrodynamic 
impact. It also presents comparisons of numerical solutions of these 
equations with theoretical fixed-trim loads, moments, and motions under 
identical initial conditions, and a comparison is made of experimental 
free-trim and fixed- trim data. 
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SYMBOLS 

hydrodynamic aspect ratio (tan p/tan t) 

distance between center of gravity and step, parallel to keel 
distance between center of gravity and step, normal to keel 
hydrodynamic force 

moment of inertia about center of gravity (pivot point) 
constant 

wetted keel length 

moment about center of gravity (bow- up moment is positive) 

moment about step 

virtual mass of float 

mass of fluid in a flow plane 

acceleration of center of gravity, parallel to keel 
distance of flow plane forward of step, parallel to keel 
time after contact 
velocity 

horizontal acceleration of center of gravity 
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x horizontal displacement of step from contact point 

• • 

Y vertical acceleration of center of gravity 

y draft at step 

• * 

Z acceleration of center of gravity, normal to keel ' 

z penetration of float into a flow plane 

0 angle of dead rise 


7 

P 


flight-path angle 



mass density of water 


r trim, radians, except as noted 

F(p) dead-rise function 

cp(A) aspect-ratio correction factor 

Subscripts: 

o condition at time of water contact 

n direction normal to keel 

p direction parallel to keel 

Where units are not specified any consistent system of units may he 


used. 


Dots over x, y, z, and r indicate differentiation with respect 
to time. 


METHOD OF AMLISIS 


Physical concepts .- The physical, concepts on which this analysis is 
based are discussed in detail in references 1 and 3> in which the solu- 
tion is restricted' to fixed positive trim. Briefly, the theory discussed 
in these references was based on the concept that the primary flow about 
an immersing seaplane float or hull occurs in transverse planes which 
were considered fixed in space and oriented perpendicular to the keel. 
Figure 1 is a sketch of a prismatic float immersing at positive trim. 
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A given flaw plane is indicated in figure 1(a) at a distance s forward 
of the step and the penetration of the float into this flow plane 1 b 
represented by ’ z. The motion of the fluid in each plane was treated 
as’ a two-dimensional phenomenon independent of the other flow planes. 

The total force on the float is obtained by integrating the reactions 
of the fluid in the individual flow planes in contact with the float. 

This force is reduced by the application of an end-flow (aspect-ratio) 
correction to account for the losses which exist in three-dimensional 
flow. The effects of buoyancy and viscosity are neglected, since in an 
impact they are normally smal 1 in comparison with the inertia forces. 

A lift force equal to the weight of the seaplane is also assumed to act 
throughout the impact. 

These concepts are followed in the present analysis except for the 
modifications necessary to introduce freedom in trim. In. order to 
introduce freedom in trim, the flow planes are assumed to rotate about 
instantaneous centers of rotation at their intersection with the float 
keel. Thus, the planes are permitted to maintain their orientation 
normal to the keel and at the same time to have no translation in space 
along the line of the keel. 

Equations of motion .- In this analysis both the float being con- 
sidered and the forces applied to it are symmetrical about a vertical 
plane through the keel so that the float motions can be resolved as a 
case of plane motion. The equations of motion for the float (see 
fig. 1(a)) can then be written as 
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• * • • 

In equations (1) and (2) the quantities S and Z are the com- 
ponents of the center-of-gravity acceleration which are, respectively, 
parallel and perpendicular to the keel. These accelerations can be 
expressed in terms of the horizontal and vertical conponents of the step 
acceleration and the center-of-gravity accelerations relative to the 
step in the following mariner: 


( 4 ) 


• • •• • • • p • • 

S = x cos t - y sin t - a-r^ - bT 
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and 


•••• •• • Q •• 

Z = y cos t + x sin t + hr* - aT 


(5) 


Since the viscous forces of the fluid are considered small enough 
to he neglected in this analysis, there is no force on the float in the 
direction parallel to the heel and hence Fp is equal to zero. There- 
fore, from equations (l) and (4) the following expression can he obtained: 


0 = -mj> (x cos t - y sin t - &t 2 - bir) 


( 6 ) 


This equation can then he solved to obtain an expression for x: 


.. y sin t a-r bT 

x = + + 

COS T COS T COS T 


(7) 


Now consider equation (2), the normal-force equation. The normal 
force on the float is the sum of the forces exerted by the individual 
flow planes reduced by an end-flow- loss correction factor. The typical 
flow plane illustrated in figure 1 (a) is at a distance b forward of 
the step and has a width ds parallel to the keel. The float has made 
a penetration z into the plane. From the analysis of reference 3> 
the virtual mass of fluid in this flow plane is 

m* = Kz 2 ds ( 8 ) 


where 


K = |p*F O) (9) 

In accordance with the assumptions, the reaction of the fluid in 
the plane is the time derivative of the momentum so that 


dF = — (i%z) 

dt 


— (Kz 2 ds z) = K(z 2 z + 2 zz 2 )ds 
dt 


( 10 ) 
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The normal force acting on the float can now he obtained by inte- 
grating over the wetted length and multiplying by the aspect-ratio cor- 
rection factor 


P 


n 



dF = Kcp(A) 



(z^z + 2zz^)ds 


( 11 ) 


where F n can be seen to be a function of z, z, and z. These quan- 
tities may be established as follows. 

It can be seen from figure 1 (a) that the penetration z of the 
float into a given flow plane can be written as a function of the vari- 
ables y, s , and t, each of which is a function of time. Thus 


z = 


r s sin t 
cos t cos T 


( 12 ) 


Differentiation of z with respect to time gives 


z = 


s sin r yT sin r bt 
+ _ 


COS T COS T COS^T 


COS 2 1 


( 13 ) 


and a further differentiation gives 


.. y s sin t y*r sin t st yj^ 2yr 2 sin 2 T 2sT 2 sin t 
z + - — : — + + + 


COS T COS T COS^T COS^T COS T COs3t COS 


3t 


_ .2 

2yT sin r 2 st 


2 ? 
COS T COB c T 


( 1*0 
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The derivatives s and s* are now expressed in terms of x, y, 
and t. In accordance with the assumptions, the velocity s is con- 
stant over the wetted length and is simply equal to the negative of the 
float tangential velocity Vp j thus, from figure 1(h), 

s = -Vp = y sin t - x cos t ( 15 ) 

The derivative of this expression is 

. . ds . *♦ • • •*, 

s = — = y sin t - x cos t + yr cos r + xt sin t 

dt 


Introducing equation ( 7 ) gives 


s = yT cos 


t + xr sin t - af^ - b*T 


( 16 ) 


Equations (12), (13) , and (l4), with s and s* replaced hy equa- 
tions ( 15 ) and (l6), respectively, are now introduced into equation (ll). 
Integration is then performed to obtain the normal force in terms of 
motions of the step. Substituting the resulting expression for normal 
force into equation (2), where *Z* in this equation is replaced by equa- 
tion (5), gives the following equation for motion of the center of gravity 
normal to the keel: 


Kcp(A) 


12 sin^T cos^t 


A{v (4 sin 3 T cos t - sin t cos t) + T 2 (l6 sin^T - 


6 sin 2 T + 2 )j + y^j^k sin 2 T cos 2 t) + V(4b sin^T cos 2 t) 


T^(ka. sin^T cos 2 t) + xt (20 sin^T cob 2 t) + yr(l 2 sin r cos^t - 


20 sin t cos5t)J + y 2 jl 2 sin 2 T cos^T(y cos t + x sin t)^~] 


>■ = 


/ y bT sin r ,p .. ai^sin 
— mpl + + br - ar + 


:os t cos T 


cos 




(17) 
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In order to solve the third equation of motion, the total moment 
exerted "by the individual flow planes about the center of gravity must 
he determ in ed. To make an approximate correction for end-flow losses 
this moment is multiplied by the aspect-ratio factor cp(A) , where the 
assumption is made that this over-all correction factor can be applied 
to the moment as well as to the normal force. Equation ( 3 ) may thus be 
written as 


M 


eg 


= 9(A) J (s - a)dF = I cg ‘i 


( 18 ) 


which, with the use of equation (ll), becomes 


oy/sin t 

M cg = Kcp(A) J (z 2 z‘ + 2zz 2 )s ds - aP n ( 19 ) 


Proceeding as with the normal-force equation gives the following equation 
for angular motion about the center of gravity: 

— "T~ c~i y ^[^(5 sin^T cos t - 2 sin t cos t) + t 2 (25 sin^ - 

60 sin^T cos^r [_ — 


19 sin 2 T + 6)J + y^|y(5 sin 2 


t cos^t) + T(5a sin 2 T cos t - 


20a sin^T cos t + 5b sin^T cos 2 t) + t 2 (35s- sin^T - 85a sin^T - 

10a sin t) + xt( 25 sin 2 T cos 2 t - 35 sin 2 T cos^r) + yr(l5 sin t cos^t - 

35 sin t cos5t)J + y3^r(20a sin^T cob 2 t) - *-r(20ab sin^T cos 2 r) - 

T 2 (20a 2 sin^T cos 2 t) _ xT(l00a sin^T cos 2 t) + yr(A0a sin 2 T cos^t - 

100a sin^T cos^x) + 20 sin 2 T cos^xCy cos t + x sin x)^J - 

y^j^Oa sin^T cos^T(y cos t + x sin t) !]} - w < 20 > 
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Equations {l) , (17) .> and (20) are the equations of motion for the 
float and involve the three unknowns x, y, and t. These equations 
can he solved numerically for the time histories of the three variables 
and their derivatives; see appendix A for a step-by-step iteration pro- 
cess. Involved in the solution of these equations are the dead-rise 
function (see equation ( 9 )) and the aspect-ratio correction factor (see 
equation (ll)), which are discussed in the subsequent section. 

Factors F(p) and qp(A) .- In order to perform the numerical solu- 
tions expressions must be obtained for the dead-rise function P(P) and 
the aspect-ratio correction factor cp(A) contained in the normal-force 
and moment equations. 

An expression for the dead-rise function F(p) , which gives the 
variation of the effective two-dimensional fluid mass with dead rise, 
was given by Wagner in reference 4 as 


F(P) = 


■M 


which is equivalent to 




of reference 3. 


Although this relationship has not been experimentally verified 
for intact s of floats having low dead-rise angles, it has been found 
to be in -substantial agreement with test data obtained with floats 

-J O Q 

of 2 2i to 1*0 angle of dead rise (see reference 3)- 

On the basis of experiments with vibrating plates in water (refer- 
ence 5 ) ) Pabs/fc derived the following aspect-ratio expression for approxi- 
mating the three-d i mensional virtual mass from the virtual mass computed 
on a two-dimensional basis: 


<p(A) = 1 - — 
2A 


where A is the aspect ratio of the equivalent vibrating plate. If 
this reduction in the virtual mass is assumed to be determined by the 
shape of the intersected area in the plane of the water surface, then 
the application of Pabst’s data to V-bottom floats results in the 
expression 

tan t 


<P(A) = 1 - 


2 tan p 
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THEORETICAL RESUETS 


Reduction to fixed trim .- The free-trim equations were derived 
from the same basic assumptions as the fixed-trim equations; they must, 
therefore, he reducible to the fixed-trim equations. This reduction is 
accomplished in appendix B by considering the trim to be a constant 
instead of a variable, so that the terms containing t and t drop 
out. It is then demonstrated that equations (17) and (20) reduce to a 
form identical with the equations for normal force and step moment in 
references 1 and 6 , respectively. Equation (7) is also shown to reduce 
to a simple fixed-trim relation. 

Comparison of free- and fixed- trim theoretical solutions .- In order 
to indicate the effect of freedom in trim both free- and fixed-trim 
theoretical solutions were made for a .typical hull for three different 
sets of landing conditions. The dimensions and inertia values of the 
full-scale float considered are given in table I. The theoretical 
initial conditions for the three impact cases treated are listed in 
table II. (The experimental values also listed in table II are discussed 
in a subsequent section. ) 

Figure 2 presents time-history comparisons of solutions of the 
free-trim equations with fixed-trim solutions for the three cases of 
theoretical initial conditions given in table II. The quantities 
plotted are the vertical and angular displacements, velocities, and 
accelerations . 


The free-trim solutions were obtained by the method described in 
appendix A and the fixed-trim solutions were obtained by the method of 
reference 7- The so-called t curves for fixed trim shown in the plots 
were obtained by using equation ( 3 ) hut with the fixed- trim total moment 
from equation (B7) instead of the free-trim moment. These r curves 
can be interpreted as a measure of the applied moment since, from equa- 
M 

tion ( 3 ), r = — £ and I C g is the constant value for the free-trim 
I cg 


float. 


EXPERIMENTAL RESULTS 


To demonstrate the effects of freedom in trim experimentally, both 
free- and fixed-trim tests of a float model were made in the Langley 
impact basin. The measured motions are compared in this section. 
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No direct comparison "between the experimental and theoretical 
results can he made because of differences between the float used in 
the tests and that assumed in the theory. The float used in the tests 
had a pointed step, flared chines, and an afterbody; whereas the theory 
applies to a float forebody having a transverse step and an infinitely 
wide V-bottom. Since these float tests represent the only free-trim 
impact-basin data now available, however, the results are included 
herein for the qualitative comparisons that can be made. 

Apparatus and instrumentation .- The Langley impact basin and testing 
procedure are described in reference 8. Pertinent dimensions of the 
float model used in the experiments are given in table I. This model, 

except for the differences outlined previously, is a scale model of 
the float considered in the theoretical solutions. 

In the free-trim tests the float model was attached to the dropping 
linkage of the launching carriage at two main pivot points; these sup- 
ports were on a transverse line which passed through the point corre- 
sponding to the center of gravity of the . airplane from which the float 
was patterned. For the fixed-trim tests a third support point was 
located about 20 inches aft of the main supports on the longitudinal 
center line of the float. 

For the free -trim tests the model was supported only at the two 
main points and was held at fixed trim until just prior to contact by 
means of a locking mechanism. After contact the model was free to 
rotate in pitch about the main supports over a trim range of -6.5° 
to 22. 5°- Beyond those limits angular displacement of the float was 
restrained by two shock struts which were coupled to the float by means 
of telescoping tubes, one 60 inches forward and one 60 inches aft of 
the main pivots. The buffer action of the shock struts extended the 
trim range about 5° in each direction before a stop was reached. 

Two strain-gage accelerometers of the same type of construction 
were electrically connected to obtain angular acceleration directly. 

These accelerometers were located on a longitudinal line passing through 
the axis of rotation and at a distance of 6 feet forward and 6 feet aft 
of this axis. A control-position transmitter was adapted to the equip- 
ment to measure angular displacement. 

A standard MCA three-component accelerometer was used to obtain 
the vertical component of acceleration of the float. It had a natural 
frequency of 21 cycles per second and a critical damping of 0.8. Other- 
wise, the standard instrumentation as described in reference 8 was used. 

Comparison of free- and fixed- trim experimental results .- Figure 3 
presents the comparison of free-trim and fixed-trim experimental results 
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for cases I and II listed in table II. Again vertical and angular 
displacements , velocities , and accelerations are compared; the values 
shown apply to full-scalq conditions. 

• It can he noted in table II that the initial vertical velocity of 
the fixed-trim run of case I was higher than that of the corresponding 
free-trim run. The time histories for the fixed-trim run were there- 
fore adjusted as follows to give more comparable results in figure 3(a). 
Since the initial flight-path angles of the two runs were not greatly 
different and the initial tr im s were the same, the velocity time history 
of the fixed- trim run was scaled by the ratio of the free-trim initial 
velocity to the fixed-trim initial velocity, the acceleration time 
history was scaled by the square of this ratio, and the time was scaled 
by the reciprocal of the ratio. No corresponding adjustments were made 
for case II. 


DISCUSSION 


Initial conditions .- Table II presents the initial conditions for 
both the free- and fixed-trim theoretical solutions and also for the 

free- and fixed-trim test runs. The test runs were made with a - - scale 

2 

model but the initial conditions are given in terms of full-scale values. 

The initial conditions for the theoretical solutions of cases I 
and II were chosen to correspond to the free-trim experimental runs of 
those cases to make the results as conparable as possible. The experi- 
mental fixed-trim runs were chosen, from among those available, to 
correspond as closely as possible to the respective free-trim runs. 

The initial conditions for case HI represent a more severe impact 
such as might occur in a second or third contact during a landing run 
or in a contact against the flank of an advancing wave. 

Case I .- The first set of conditions treated (case I) represents 
an impact at moderate initial values of flight-path angle and trim. In 
both the theoretical (fig. 2(a)) and the experimental (fig. 3(a)) results, 
the vertical motions, trim, and angular velocity for the free-trim con- 
dition differ only slightly from the corresponding motions in the fixed- 
trim condition. The only angular acceleration or moment of any size is 
exhibited by the experimental free-trim run (t ~ 0 . 15 ), but this moment 
is associated with the immersion of the flared chines. Thus, for the. 
initial conditions of case I no significant difference appears to exist 
between fixed- and free-trim impacts. 
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Case II .- Case II has approximately the same initial flight-path 
angle as case I hut has a .higher initial trim. The vertical motions 
for the free-trim condition are again observed to approximate those 
for the fixed-trim condition in both the theoretical (fig. 2(b)) and 
experimental (fig. 3(b)) curves. This impact, however, resulted in an 
appreciable change of attitude and a moderate nose-down moment. 

The initial trends of the theoretical and experimental angular 
motions are similar. However, the immersion of the flared chines prob- 
ably caused large loads on the float forward of the center of rotation 
and therefore a large, nose-up moment as indicated by the positive peak 
of the experimental r curve. These loads would result in large 
deviations of the subsequent motions of the test float from those of 
the float assumed in the theory. 

Case III .- Case, III represents a high- flight-path-angle, low- trim 
condition such as might occur in an impact on .the flank of a wave. The 
curves represent the "effective" values referred to the inclined plane 
of the water surface as in reference 9 and are analyzed as an extreme 
condition of smooth-water impact. 

In this case, the deviations between the fixed- trim and free-trim 
results are much greater than those previously described; freedom in 
trim in fact reduces appreciably the maximum values of y and r. The 
large values of negative vertical velocity of the free-trim solution 
compared with the fixed-trim solution late in the impact appear to be 
consistent with the large values of trim attained. 

Probably the most serious effect observed for the free-trim solu- 
tion is the high trim and large positive value of angular velocity at 
the end of the impact. This could result in a stalled condition between 
impacts or could lead to extreme initial conditions for a subsequent 
impact. However, these conditions should be somewhat restricted by the 
presence of an afterbody and by aerodynamic moments not taken into 
account by the present theory. 

No impact-basin data have been obtained for conditions as severe 
as those in this case. An indication of the miotions to be expected, 
however, can be obtained from an impact of a four-engine flying boat, 
the data for which have not been published. The initial conditions for 
this impact were substantially the same as those of case III and the 
resulting motions were also very similar. Of particular interest in 
this seaplane impact is the increase in trim from a small initial value 
(about 3°) to about 15 ° with positive angular velocity at the time of 
exit from the water. The acceleration records contained large oscilla- 
tions, due probably to structural vibrations, which rendered them 
unsuitable for direct comparison with the theory; however, the peak 
values of the faired curves for both y and t agree roughly with the 



MCA TN 2698 


ik 


corresponding values for the free-trim solution of case III. Thus, 
although direct comparison is not practicable, this flying-boat impact 
indicates that the initial conditions for case III are within the 
practical range and that the theoretical results which are obtained 
for this case are at least consistent with flight experience. 

General discussion .- In general, for cases I and II, the theoreti- 
cal free-trim vertical motions 'approximate the fixed-trim vertical 
motions to about the same degree as was found experimentally. Also, 
the trends of the theoretical angular motions are similar to those of 
the experimental data except where deviation is to be expected because 
of differences in body shape. Thus, for impacts having moderate initial 
conditions the theoretical and experimental free-trim results appear 
to be in agreement. 

For the more extreme ini tial conditions of case III, there is, of 
course, more deviation of all variable's from the fixed-trim solution. 

The deviations which occur are, however, qonsistent among themselves 
and with general experience. Also, the reductions of maximum applied 
load and moment which occur are to be expected. 


COKCHJDING REMARKS ' 


Equations of motion are derived and presented for free-trim hydro- 
dynamic impact of a V -bottom, transverse step, prismatic float forebody. 
The equations are an extension of previously published fixed-trim hydro- 
dynamic impact theory and are chiefly based on the same concepts and 
assumptions. A method of solution of the equations is also presented 
which gives the results as time histories of vertical, horizontal, and 
angular displacement, velocity, end acceleration. 

The free-trim equations are shown to be reducible to the fixed- 
trim case. Moreover, comparisons of specific solutions of the free- 
trim equations with fixed-trim solutions for the same initial conditions 
and also with some experimental data have shown that the computed 
results are reasonable. The free-trim solutions and experimental data 
exhibit similar trends and deviations from the fixed-trim case, and 
changes of attitude and applied moment are consistent with previous 
experience. 

Although these facts do not completely validate the equations, 
they at least indicate that reasonable and consistent results can be 
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obtained by their use. Further experimental investigation should 
provide a better evaluation of the method. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., February 13, 1952 
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APPENDIX A 
METHOD OF SOLUTION 


Before the numerical solution can he performed, equations (17) 
and (20) must he transformed into forms better suited to the method of 
solution used. Equation (7) is already in the form required for 
solution. 

Equation (17) can he rewritten as 


and equation (20) as 


A*t + By = C 


• • • • 

Dt + Ey = F 


(Al) 


(A2) 


where 


A = 


_ K<p(A) [A 4 sin3 T cos t - sin r cos t) + y3(4b sin^T cos^^J 

■12 sin^r cos^t 


/ h sin t \ 

m fl a) 

\ COS T / 


(A3) 


B = 


_ Kcp(A) [y3(4 sin^T cos^t)] ^ 


12 sin^T cos^t 


COS T 


(A4) 
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K<P(A) 


3 ^ ( 16 sinV - 6 sin 2 ! - + 2)] 


12 sin^T cos^r 

sin 3 r cos^r) + if (20 si n k r cos 2 -!-) + if (12 sin t cos 3 t 


20 sin t cos 


^ + y 2 jl2 sin 2 T cosS^y cos t + x sin T)^jj* “ 


( aT 2 sin t . 

+ br, 

. COS T / 


(A5) 


D = 


— — [y 5 (5 sin 3 T cos t - 2 sin t cob t) + 

T COS-^T — 


Kcp(A) 

[ 

60 sin^r 

cos 5 t *- 


y^(5a sin 2 T cos t - 20a sin^T cos t + 5b sin3-r C os 2 t) _ 

y3(20ab sin^T cos 2 t) I - I 

I C S 


(A6) 


E 


Kcp(A) 


60 sin^ 


sin 2 T cos 2 t) - y3 ( 20a sin 3 r cos 2 t)~] 

T COS^T — I 


(A7) 
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K<p(A) 


60 sin^r cos^t 


- y 5^(25 sin^T - 19 sin 2 T + 6 )J + 


i 2 ( 35 a sin^T _ 85a sin^T - 10 a sin t) + xr(25 sln^T cos 2 t - 


35 sin 2 T cos^t) + yr (15 sin r cos^t - 35 sin r cos^t)J - 


n 


3 T 2 ( 20 a 2 sin^T cos 2 r) + xV(l00a sin^r cos 2 t) - yr(itOa 8 in 2 T cos^t - 


I 


100 a 


sin^r cos^t) - 20 sin 2 T cos^r(y cos t + x sin r)^j - 


y 2 [60a sin^T cos^i^y cos t + x sin r) 2 


(A8) 


Equations (Al) and (A 2 ) can lie solved simultaneously to give 


.. AF - DC 

y = 

AE - DB 


(A9) 


and 


•• EC - BF 

T = 

EA - BD 


(A 10 ) 


which, together with equation (7), 


• • • 2 • • 

. . y sin t ar hT 

x = + + 

COS T COS T COS T 


are the equations to he used in the numerical solution. 
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The numerical solution can "be performed hy extrapolation of the 
plots of y and t against time by any desired function to obtain an 
approximate value at the end of the next time interval. The other quan- 
tities (x, x, x } y , y } Tj and r) are then confuted and the process 
iterated. By this method, time-history plots of each of the nine vari- 
ables are obtained. The present solutions were made, however, by the 

Kutta ^ method (reference 10) which is more readily adaptable to the • 

automatic computing machine used. 
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APPENDIX B 

REDUCTION TO FIXED-TRIM CASE 


One test of the validity of the free-trim, equations is their 
conqmtibility with the special, and well-established, case of fixed- 
trim intact. If, in equation (17), the trim is considered to remain 
constant throughout the impact, all terms containing derivatives of 
the trim are eliminated and the equation becomes 


Kcp(A) 


12 sin 3 T COS^T 


%"3y sin^T cos^t + 12y^sin^T cos^T(y cos t + x sin t)^ = 


-m-p 


• • 

y 

COS T 


(Bl) 


which reduces to 


m f y Kcp(A)y3y Kcp(A)y 2 (y cos t + x sin t ) 2 

= 5— + 

cob t 3 sin t cos^t sin t cos t 


(B2) 


However, 


= z = 


cos T 


dt 


(B3) 


and 


Y n = y cos t + i sin t ' (B4) 


so that, by introducing equations (2), (B3), and (b 4) into equation (B2), 

dV 

Kcp(A)y 3 ^ Kcp(A)y 2 V n 2 

F n = — + (B5) 

3 sin t cosS- sin t cos t 


which corresponds to equation (22) of reference 1. 
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Similarly, the right-hand side of equation (20) is equal to the 
moment about the center of gravity as in equation (3) and, if the trim 
derivatives are set equal to zero in the left-hand side', the equation 
becomes 

— — f^y^y sin^T cos^t - 20ay3y sin^r cos^t + 

60 sin^ cos^t *— 


20y3g±n^T cos^T(y cos r + x sin t)^ - 


60ay^sin3r cos\-(y 


cos r + x sin 




■which reduces to 


K<p(A)y^y Kcp(A)y3(y cos t + x sin r) 2 Kcp(A)ay3y 


M—rr — " ■ ■ — - + 

Cg p O 

12 sin t cos j t 


3 sin^r cos r 


3 sin t cob-^t 


KcpCAjay^y cos r + x sin t)‘ 


sin t cos t 


After equations (B3) and (BA) are substituted into equation (B7), the 
moment becomes 


Kcp(A) y^y -3P/- ♦ ^2 

M-o. — + y^cos Tly cos t + x sin t) 

c b 2 3 A 

3 sin t cos j t1_7 . — 


_ U.V 

K<p(A)y3 Kcp^y 2 ^ 


3 sin r cos 2 t s ^- ri T cos T 


and, by introducing equation (B5) and equations ( 9 ) and (ll) of refer- 
ence 3, equation (B8) becomes 



22 


MCA TN 2698 


Meg + ^n 


Kcp(A) 

A* • 

y t 

3 sin^r cos^t 

_k 


y3(y + K ± 


cos t) 


2 


(B 9) 


-where K-j_ is K In reference 3- 

From equations (3^-) and (2 6 ) of reference 6 the moment about the 
step is 

' M s = M cg + aF n (BIO) 


Substituting equation (BIO) into equation (B9) gives the moment about 
the step 


M s = 


Kcp(A) 

y^y 

3 sin^T cos^t 

Jj- 


+ y 3 (y + Ki cos t)' 


(Bll) 


which is identical with equation (l 6 ) of reference 6 . 

Equation (7) can also be reduced in this maimer to give' 

x = y tan t (B12) 

In the fixed-trim condition the step motions and center-of-gravity 
motions are the same. Equation (B12) is, therefore, equivalent to 

' X = Y tan t (B 13 ) 


Equation (B13) can readily be deduced from figure 1(c). 
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TABLE I 

1 - i 

PHYSICAL CHARACTERISTICS OF FULL-SCALE HULL AND =r - SCALE ; 

2 . \ 

) 

EXPERIMENTAL MODEL 


Experimental 



Full-scale 

model 


hull- 

- scale^ 

Beam, ft 

. 6.33 

3.17 

Forebody length, from centroid of 'step, ft . 

. 19.11 

9.56 

Over -all length, ft 

. 35.67 

17.83 

a (from centroid of step), ft 

2.28 

l.l4 

h (normal to keel), ’ft 

6 . 06 

3.03 

Weight, lb . 

lk, 000 

1,750 

mf, lb-sec^/ft 

434.78 

5^.35 

I C g, lb-sec^-ft 

19,43,0 

607 

P, deg 

20 

20 



1 



TABLE XI 


INITIAL CONDITIONS FOR TIME HISTORIES 


paae I 


Case II 


y 0 > ft 

y 0 > ft/aec 
*y 0 , ft/sec/sec 


Xqj ft/sec/sec 

f radian b 
T °^\ deg 

t 0 , radiana/sec 

t q , radians/eec/sec 

7 0 > de B 


Theoretical 
free and 
fixed- trim I Free 


Experimental 

(a) 


Fixed 

trim 


120.4 125.7 


0.122 0.122 

7.0 7-0 


Theoretical 
free and 
fixed trim Free 
trim 



Experimental 

(a) 


1 £ 3.7 


0.209 

12.0 


Fixed 

trim 


125.7 


Case III 


Theoretical 
free and 
fixed trim 



, 4w naca in 2698 






























ft/sec 2 7, ft/sec 






Tine after contact, t, sec 





, ft/sec 2 y, ft/see 
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